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Abstract 

We investigate the instabilities of low winding number electroweak strings us- 
ing standard numerical techniques of linear algebra. For strings of unit winding 
we are able to confirm and extend existing calculations of the unstable region in 
the (m,H /mw, sin 2 Oyy) plane. For strings of higher winding number we map the 
unstable regions for the various decay modes. 
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It has been known for some time that certain field theory models with spontaneous 
symmetry breaking contain toplogical defects, such as domain walls, strings and mon- 
poles. These defects are stable because their field configurations wrap around the vacuum 
manifold, and to unwind them requires lifting the Higgs field to its vacuum value in the 
unbroken phase, throughout an infinite volume of space. There is therefore an infinite 
energy barrier separating these defect solutions from the vacuum. 

Recently, defect solutions have been found in theories where they are not topologically 
stable jl], 0, 0]. These solutions may be a maximum or a local minimum of the energy, 
depending upon the dynamics of the model in question. Such 'embedded defects' have 
been found in the electroweak model: SU{2)l x U(1)y — > U(l) em . This contains two 
embedded U(l) — > I symmetry breakings which can both give rise to Nielsen-Olesen 
vortex solutions ||. These solutions are called r- and Z-strings depending upon how 
one chooses the embedding, r-strings are expected to be unstable for all values of the 
parameters M and so we will only consider Z-strings. 

The stability of Z-strings for a range of parameters was addressed in James et al 
H where it was shown that for physical values of the parameters the N = 1 string is 
unstable. The stability analysis did not extend down to small values of the Higgs self 
coupling because of numerical problems caused by the widely different scales of the two 
cores of the string. The method adopted for fixing the gauge, locating the unstable mode 
and computing the eigenvalues was not wholly transparent either. 

In Perkins || and Achucarro et al |7| instabilities due to ^-condensation were in- 
vestigated, and in Achuccaro et al it was shown that a iy-condensed string of winding 
number N was equivalent to a string of winding number N — 1. In addition to this decay 
mode Achuccaro et al also showed that there is another distinct physical decay mode 
for large N strings, where the upper component of the Higgs field acquires a non-zero 
value at the string core. This decay mode should result in the complete unwinding of the 
string. 

MacDowell and Tornkvist M showed that for some, unspecified, regions of the pa- 
rameter space a string of winding number N is unstable to the formation of VT-fields 
with angular momentum m such that — 2N < m < 0. 

In this letter we seek to complete the stability line of James et al, investigate the two 
decay modes of Achucarro et al and investigate the 2N — 1 decay modes of MacDowell and 
Tornkvist. We are careful to fix the gauge at the outset, and we use more usual eigenvalue 
methods than were adopted in James et al. We do not, however, follow Achucarro et al 
in tracking the full non-linear evolution of the fields. 

We consider the bosonic sector of the Weinberg-Salam model with the lagrangian 

C = -\b^ - \G%G a ^ + {D^)\D^<j>) - A(0V - n 2 /2f, 

where 

D, = {d (t -iga a W^/2-ig'B lt /2)<f> 
Bau = duB v — d u Bu 
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a a are the Pauli spin matrices, are the SU(2) gauge fields, are the U{1) gauge 
fields and cf) is the Higgs doublet. The SU(2)l x U(1)y symmetry of this lagrangian 
is broken to the U(l) em symmetry of electromagnetism by the Higgs field acquiring a 
vacuum expectation value <ft = (0,i]/^/2). With the standard field basis 

= cos 9 W W 3 - sin dwB^ = sin 9 W W 3 + cos OwB^ 

K = - iW D w » = ^« + iW D 

the boson masses generated by the Higgs mechanism are 

1 1 



M w = -g V , M z = -g zV , M H = ^2\^. 



The Weinberg angle 6w is given by taxiOw = g' /g, and g z = \jg 2 + g' 2 . 

Consider time independent field configurations with (f> u = = W~ = = and 
Zq = 0, then the energy functional is given by 

E = Jd 3 x {^ZijZij + |(Vi + '-fzM 2 + A(|0 d | 2 - r, 2 /2) 2 ^ , 



where 



Zij — V 1 %Zj — WjZi. 



This is the same as for the Abelian-Higgs model and so there is a Nielsen-Olesen string 
configuration given by 

n 2N 
<t>d=^f(r)e me , Z r = 0, Z e = a(r), 

which extremises the energy. The functions f(r) and a(r) are identical to those found by 
Nielsen and Olesen [[J. 

Consider a general perturbation to the solution of the form 

K = e56 u e- iujt d>* ; = e56*e iujt 



u 



W+ = e5W+e~ iuJt W~ = ebW~e~ iuit 

Z^ = Z' + eSZpe-™ A^ = e5A fM e^ t . 



This gives corrections to the action 

S = S(4>% ZD + e 2l -J d A x5 X ] V5 X + 0(e 3 
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where T> is the perturbation operator and 



5x ] = (Sfae-**, Sfce™, SW+e^, 5W~e~^\ Sfae-™, 5Z,e~^\ SA^). 

The equations of motion for the perturbations T>5x = 0, could in principle be coupled 
equations in twenty fields. For the string background however, they separate into four 
sets of coupled equations. The first two are of the form 

which, since the string solution is invariant under charge conjugation, are just charge 
conjugates, so we need only consider the first of these. The third is of the form 



',/ I - o, 



\sz, 



which contains only the fields of the Nielsen-Olesen vortex, and was addressed in ||10|| . 
The fourth equation is 

V,5A^ = (rT<9 2 - rd*)6Ap = 0, 

which is just the usual wave equation of electromagnetism. The last two equations do 
not contain any modes which can decrease the winding number, and so we must look at 



the solutions of (1) above for decay modes. As mentioned in [IIJ these equations contain 
linear derivative terms and gauge degrees of freedom, both of which can be removed by 
choosing the background gauge conditions 

Fl (W+) = d»W+-ig cos 9 W Z2W+-^?0 U =O (2) 

F 2 (W~) = d^W; + ig cos 9 W Z^W~ + ^=^1 = (3) 

F 3 (Z„) = ^-^(0^-C0 d ) = O (4) 
F 4 (A,) = 0M„ = O (5) 

which are imposed by adding the gauge fixing terms 

1 4 2 

i=i 

to the lagrangian. This enables us to separate out the linear time derivatives and to set 
up gauge- fixed eigenvalue equations, which for equation (1) above are 

(«) 
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and 



(-V 2 + 2ig cos 9 w ZtV k + g 2 cos 2 9 w Z^ k + ^| 2 ) ( ^ ) = ^ ( SW+ ) ' (?) 
where i = 1,2, and the gauge fixed perturbation operator is given by 



r <;r ( Dl Xi 

X* Do 



where 



2 2 

D 1 = -V 2 +tg z cos29 w Z c k V k + ^cos 2 29 w Z k "Z c k + 2\(\ ( p d \ 2 -r ] 2 /2) + Ucf> d 



a 2 

D 2ji = M-V 2 + 2ig cos0 w Z£V* + g 2 cos 2 9 W Z c k Z c k + ^-\(f) d \ 2 ) + 2ig cos 9 W Z 



X t = tg^Vt + ^Z^ 
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x; = -igV2(y 3 - l -fz<w. 

The 5Wq and 5W+ perturbations decouple because the background string solution is 
independent of t and z. 

The gauge condition used means that the gauge fixing terms must be accompanied 
by Fadeev-Popov ghost terms. This results in an accompanying eigenvalue equation for 
the ghost field A + 

a 2 

(-V 2 + 2ig cos 9 w Z c k V k + g 2 cos 2 9 w Z c k Z c k + ±-|<^| 2 )A + = u 2 A + . (8) 

This is the same equation as (7) and so A + , 5Wq~ and SW^ all have the same eigenvalue 
spectra. Of the total of 5 eigenmodes in (6), (7) and (8), only three should be physical, 
corresponding to the three spin states of the massive W boson. This is ensured by the 
ghosts canceling one linear combination of 5W$ and and one of the eigenmodes of 

M GF . 

It may seem that the gauge choice (2)- (5) introduces unnecesary complications and 
a difficulty in identifying the physical modes. However, we are looking for decay modes 
where u 2 < and the eigenvalues of (8) are positive, so there is no difficulty in identifying 
the field configurations of the decay modes. It is also worth pointing out that (2)- (5) 
results in a considerable simplication in the resulting eigenvalue equations (as opposed 
to the temporal gauge for instance). 

We expand the scalar and gauge fields in total angular momentum states. For the 
gauge fields they are 

= - and W+ = ^(W;+ l -W+). 



The total angular momentum operator for the gauge fields is J z = L z + S z where 



with (S z W + )i = +W^, and (S z W + ) l = -Wf . So the suffices (t) and (|) identify gauge 
fields with spin up and spin down respectively. The background string solution in total 
angular momentum states is 

and the perturbations are 

5(j) u = ve™' 9 , SW^ + = ~iw m e 

with N + m = m'. It is useful to rescale 



,i(m— 1) 



= E^-m e * (m+1)9 ' 



2^2 



(/» = — 6 Z s = —Z s W + = —W + W~ = —W~ r 

to give dimensionless variables. The eigenvalues are now in units g 2 r] 2 /8 and the magnetic 
field strength is in units g z rf /4 = m 2 z /g z . Substituting in the above gives the eigenvalue 
equations 



/ £>i A B \ I s 



A D 2 
B D 3 



W m 
V W-m 



= uo 2 



m 1 



( 3 



(9) 



where 

D 1 = 

D 2 = 

D 3 = 

A = 

B = 
and = 8\/g 2 z . 



-V 2 + 



-v 2 + 

v p ^ 



, , v „,' + aWM 2 +w2 _ 1) + 2/W()iy 
_ + «m-l)+2aiW^ + cQs2 % + 4 cos2 iV da 

((m + 1) + 2aiVcos 2 M 2 , ,2 2 Q A 2 Q N da 

— — + 2f 2 cos 2 9 W - 4 cos 2 9 W - 

p z p dp 



2(v p /-^(l-a)) 
-2 (v„/ + ^(l-a)) 
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If we resolve s m /,w m ,w*_ m into real and imaginary parts, the complex eigenvalue 
problem seperates into two eigenvalue problems with explicitly real fields 








( 










w r 

m 


W -m J 




\ 


w r 

—m 



and 






( S m' 




( S m' \ 










) 






Win 




V - W -m J 







where D 1 ,A and B are given above. The profiles for the Nielsen- Olesen vortex were 
solved for with a relaxation method on the energy functional, and then substituted into 
the eigenvalue equations above. This eigenvalue problem was solved for a range of (3 and 
6w with the boundary conditions s m ',w m ,W- m — ► as p — ► oo by the standard matrix 
packages incorporated into matlab. 

For 1 > (3 > 0.04 a linear discretisation of p was used. For (3 < 0.04 the scales of the 
Higgs core and the gauge field core are too different to use this linear descretisation so 



we followed Yaffe [11] in using the convenient transformation 



\n{mnjmz) Vl + m^r/ 

This maps the line 0<r<ooto0<£<l. Linear discretisation of £ results in more 
points in the core region and is better suited to studying very small and very large values 
of (3 = m H /m z . With this map we were able to look at (3 down to 1 x 10~ 4 . 

The resulting stability lines are shown in the Figure. All the lines for the various 
decay modes separate the region on the left hand side where the string is unstable, from 
the region on the right hand side where the string is stable with respect to a particular 
decay mode. 

The two decay modes of Achucarro et al correspond to m' = and m = — 1. For 
the m' = mode <p u acquires a non-zero value at the core of the string, and then there 
is nothing to stop the string from completely unwinding and decaying to the vacuum. 
For the m — — 1 mode acquires a non-zero value at the core of the string, where 
its magnetic moment is aligned with the magnetic field of the string. The term giving 
rise to this TV-condensation is —4 cos 2 9w5W^~Na! / p, which is of the generic form — m.B 
for the interaction between a magnetic moment m and a magnetic field B. This term 
comes from the non-abelian nature of the field tensor. Achucarro et al showed that a 
W-condensed string of winding number N has the same energy as a string of winding 
N — 1, and so the m = — 1 mode corresponds to a unit unwinding. It is this m — — 1 
mode which is refered to as W-condensation in 0] and 0. In the background gauge, 
however, far from the string <b u can readily be indentified as the longitudinal W gauge 
field, and so all the modes we have found can be called VT-condensation. 
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For the N — l string the relation N + m = m' gives that the m' — and m — — 1 
modes are one and the same with both <p u and W\~ acquring non-zero values at the core 
of the string. The upper part of the stability line in the Figure agrees with that found 
by James et al. This is different from the stability line in || (y/]3 = 4 cos 6^) which was 
found analytically by using fairly crude approximations for the string solution and for 
the profiles of the W-fields. A large part of this difference could be that the decay mode 
of the N = 1 string has non-trivial field configurations for the M^-fields and for <p u as 
well, whereas the analysis in IJ contained no <p u field. 

For the N = 2 string as well as the two decay modes mentioned above there is a 
third decay mode for which m' = —1, m = —3. We conjecture that this decay mode 
corresponds to the N = 2 string decaying to a string of winding number N — —1. Since 
the Z-strings are non-topological and an N = — 1 string has the same energy as an N = 1 
string, such a decay is to be expected, albeit with a larger energy barrier than the decay 
to an N = 1 string. This later point is consistent with the wl — — 1 decay mode having a 
much larger region of stability in the parameter space than the decay mode to an N = 1 
string. 

Generalising this to strings of winding number N, means that we would expect there 
to be 2N — 1 decay modes corresponding to the string decaying to strings of winding 
number N — l,..., —N+l. Furthermore we would expect that the stability lines for these 
decays to be such that those for larger units of unwinding would have greater regions of 
stability, due to the increasing energy barrier for such decays. 

The stability lines for the N = 3 string, appear to be consistent with this inter- 
pretation (as are those for N—l and N = 2 strings). It is also consistent with the 
result of MacDowell and Tornkvist that strings of winding number N have regions in the 
parameter space for which they are unstable to the formation of P^-fields with angular 
momentum m such that — 2iV < m < 0. 

For the strings we studied (N = 1, 2, 3), all have a stability line close to the axis 
for the mode with m = —2N. For all the strings this line is in approximately the same 
place (sin 2 6 W ~ 0.02) and approaches the axis as the number of points in the lattice 
is increased. We believe that in the continuum limit this m = —2N mode will have zero 
eigenvalue along the line sin 2 6w = 0, and the deviation of the line from this value for 
finite lattices (100 points) is an indication of the error. 

For the complex conjugate problem with the fields <p* u and iy~,the stability lines are 
of course the same. The l W- condensation' mode is where W-T is non-zero at the core of 
the string, so the magnetic moment is still aligned with the magnetic field of the string. 

It must be noted that all the modes involve non-trivial field configurations for <j) u ,W^~ 
and W^. 

Finally, we consider the large N limit. If we assume that the U and W perturbations 
drop to zero whilst inside the core, we need only consider the small p behaviour of the 
string solution, which is / — > /nP N , a p 2 /2 as p — > 0, where fx is a constant. 
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Substituting this into the perturbation equation (9) gives 



= _ + k +i pn cos 26 w /2y + _ +2/V . cos2fa 



ft = _ v . ((m - 1) 4- , iV cos M + cog2 ^ + 4Jy cog2 g 



11' 



((m+ 1) + p 2 N cos 2 W ) 
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Ai = -V'+ vv " ' ~ y ' \ ' +2f N p 11 " cos z 6 w -4Ncos z 6 w 

A = f N Np N+1 

B = f N Np N+1 -4f N Np N -\ 

So for large N, the perturbations in cf) u and the W fields decouple. Considering terms up 
to p 2 it was shown in that the equation for <p u for m' — has a stability line given by 
1 1 — 2 sin 2 6>vf I • The Fi gure seems to indicate that the m' = line is approaching 
this limit but there is no stability region in the bottom left of the plot, which is not 
altogether suprising since iV = 3 is not large N. 

Lastly, we present an explanation for the movement of the m = — 1 stabilitiy line 
towards the line sin 2 9w = 1 for increasing winding number. Consider the equation for 
with m = — 1 for large N, with the assumption that the perturbation drops rapidly 
to zero within the core, so that we need keep terms only up to 0(p 2 ): 

' ^ - ~ + N 2 p 2 cos 4 6 W -AN cos 2 9w) w[ = uj 2 w[. 



\ dp 2 p dp 

This has a solution w\ = cexp(— N cos 2 9wp 2 /2) for which ui 2 = — 2N cos 2 9w, and so 
large iV strings are unstable to W-condensation' for all values of the parameters. This 
is consistent with the stability regions of the m — — 1 modes displayed in the Figure. 
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Figure captions 

Figure 1: The stability lines for a) the N — 1 string (solid line), b) the N = 2 string 
(dashed line) and c) the N = 3 string (dotted line). The angular momentum m of 
the modes are (reading from right to left) a) m — —1, b) m — —1, —2, —3, c) m — 
-1,-2,-3,-4,-5. 



11 



